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Abstract 

We compare different strategies aimed to prepare an ensemble with 
a given density matrix p. Preparing the ensemble of eigenstates of p 
with appropriate probabilities can be treated as 'generous' strategy: 
it provides maximal accessible information about the state. Another 
extremity is the so-called 'Scrooge' ensemble, which is mostly stingy 
to share the information. We introduce 'lazy' ensembles which require 
minimal efforts to prepare the density matrix by selecting pure states 
with respect to completely random choice. 

We consider two parties, Alice and Bob, playing a kind of game. 
Bob wishes to guess which pure state is prepared by Alice. His null hy- 
pothesis, based on the lack of any information about Alice's intention, 
is that Alice prepares any pure state with equal probability. Then, the 
average quantum state measured by Bob turns out to be p, and he has 
to make a new hypothesis about Alice's intention solely based on the 
information that the observed density matrix is p. The arising 'lazy' 
ensemble is shown to be the alternative hypothesis which minimizes 
the Type I error. 
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Introduction 



Consider two parties, Alice and Bob, playing the following game. Alice 
prepares a pure quantum state according to certain random strategy and 
then sends it to Bob. Initially Bob possesses no information about Alice's 
strategy and thus assumes that Alice performs a completely random choice 
of pure state, we refer to this statement as a null hypothesis. In this case the 
average density matrix received by Bob would be be proportional to identity. 

Measuring the received states. Bob realizes that the average quantum 
state emitted by Alice is p. However, there are infinitely many ensembles 
which average to p, and Bob still can not recover the strategy of Alice. 
Although Bob now possesses some information about Alice's intentions: if the 
received density matrix p differs from identity. Bob has to make an alternative 
hypothesis. To specify such a hypothesis, some extra principles must be taken 
into account. These principles should capture the type of Alice's behavior. 

We might assume that the strategy of Alice is to prepare eigenstates of 
p with given probabilities, but this is just an assumption that Alice is 'gen- 
erous' in providing the accessible information. Or, conversely, Alice might 
be stingy with the information and thus chooses pure states according to 
Scrooge distribution pp. 

In our game. Bob is reluctant to change his opinion and chooses among 
Alice's strategies (which average to p) the closest to his null hypothesis. By 
'closest' we mean minimizing the KuUback-Leibler [2^ distance between the 
distributions. This distance is the average likelihood ratio and is associated 
with the probability of the Type I error^. 

Another way for Bob's reasonings is to assume Alice to be lazy in efforts 
to prepare the ensemble. These efforts are quantified in terms of differential 
entropy. Remarkably, as we show in Section Q both approaches yield the 
same ensemble Q. 

1 Differential entropy and likelihood ratio 

First we have to specify a yet vague notion of 'preparation efforts' for an 
ensemble. Following [3] we formulate it in thermodynamic terms, namely, 

^To make Type I error means to accept the alternative hypothesis when the null hy- 
pothesis is still valid. An example of Type I error would be to conclude that the defendant 
is guilty, when in fact he or she is innocent. 
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we quantify these efforts by the difference between the entropy of uniform 
distribution (that is, our null hypothesis) and the entropy of the ensemble^ 
in question. The only obstacle may occur is to define this entropy, let us 
dwell on it in more detail. 

The entropy of a finite distribution {pi} is given by Shannon formula 

S{{Pi}) = ~^Pi\npi 

This expression diverges for any continuous distribution: we approximate a 
continuous distribution fi{x) by a discrete one {pi}, calculate its Shannon 
entropy, but it tends to infinity as we refine the partition. However, we 
are always interested in the difference between the entropy of the uniform 
distribution and the distribution ^{x) rather then the entropy itself. At each 
approximation step we calculate this difference, and the appropriate limit 
always exists. To show it (see, e.g. [S] for details), make a partition of the 
probability space by sets Aj having equal uniform measure. Then the 
difference between the entropies read: 

En = InN - (^-J2Pi^^P^) 

where pi = p{x)dx. The limit expression \imN->oo E^ is the differential 
entropy 

S{ii) = J /i(x) ln/i(x) dx (1) 
Remarkably, this is equal to Kullback-Leibler distance j2] 

/fi(x^ 
fi{x) In — r~^dx 

between the distribution jj,{x) and the uniform distribution fio{x) with con- 
stant density, normalize the counting measure dx on the probability space so 
that /io = 1- This distance is the average likelihood ratio, on which the choice 
of statistical hypothesis is based. Then, in order to minimize the Type I error 
we have to choose a hypothesis with the smallest average likelihood ratio. 



^We are speaking here of mixing entropy 4 of the ensemble rather than about von 
Neumann entropy of its density matrix. 
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2 'Lazy' ensembles 



The main problem reduces to the following. For given density matrix p find 
a continuous ensemble p, having minimal differential entropy ((H): 

S{p) min, j P^/i(^) dil) = p (2) 

where dip is the unitary invariant measure on pure states normalized to inte- 
grate to unity. When there is no constraints in (j2)), the answer is straightforward- 
the minimum (equal to zero) is attained on uniform distribution. To solve 
the problem with constraints, we use the Lagrange multiples method. The 
appropriate Lagrange function reads: 

C{p) = Sip) - Tr A P^piij) d^ - p^ 

where the Lagrange multiple A is a matrix since the constraints in are of 
matrix character. Substituting the expression (Q) for S{p) and making the 
derivative of C over p zero, we get 



g- TrBP,, 



m = (3) 

where B is the optimal value of the Lagrange multiple A which we derive 
from the constraint Q and the normalizing multiple 

Z{B) = [ e-T'^^P^ d^ (4) 



is the partition function for (jH)). Substituting he resulting density (jSj) to the 
expression (PJ for differential entropy we get 

S = TiBp - \nZ (5) 



3 Explicit expressions 

First evaluate the partition function Q in the eigenbasis of B. This integral 
is a special case of the calculations carried out in jH], according to which 
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Z{B) reads: 



Z(B) 



[n 



(6) 



where bk are the eigenvalues of B. If two or more of them are equal, the ap- 
propriate expression is obtained as a limit starting with unequal eigenvalues. 
To write down the expression for the eigenvalues As of the density matrix p 
via B we could evaluate the integrals 



Xs = {es\ J P^fi^tp) di)\es) 
in the eigenbasis of p. Although, like in thermodynamics, we have 



P 



dlnZ 
dB 



(7) 



which gives the explicit expression for the eigenvalues of the density matrix 



( 



n ip.-h) 



+ E 



A., 



1 



n {bs-hi) 



+ - 



E 



fc=l n (bk-bj) 



from which we see that the resulting density matrix p remains unchanged 
when we add a constant to all bk-s. That means that the matrix 'temperature' 
parameter B for the lazy ensemble is defined up to an additive constant (in 
contrast with classical thermodynamics). 

Like in [I], the expression © for the partition function can be given the 
following integral form 



ZiB) 



(n-l) 



dz 



271 i / det(S - zT) 
where the contour encloses all eigenvalues of B. 



(9) 
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So, given a lazy ensemble (jH)) with the parameter B, we have written down 
the expression (|7j) for its average density matrix. This expression is well- 
defined for any matrix B. The existence problem remains: given a density 
matrix p, is there a lazy ensemble with appropriate parameter B which aver- 
ages to p? Similar question — the existence of temperature function — arises 
in thermodynamics. The idea to solve it is the following [Sj: we consider 
the n-dimensional CDF (cumulative density function) of the measure p and 
study the asymptotics of its Laplace transform. As a result, it can be shown 
that B exists for any full-range density matrix p. 

4 Special case: spin- 1/2 particle 

In this case the state space has dimension 2. Write down the parameter B 
in the eigenbasis of the density matrix p in a suitable form: 

B = *-i+(^';a) (10) 

Then the expression (JH)) for partition function reads: 

_fe e'^ — e~^ _^ sinh /3 



2/3 /3 



(11) 



Calculating the partial derivatives according to (jH)), we get the following 
expressions for the coefficients Ai^2 of the density matrix 

where 



X,, = -±-[cotl., 3 --^J^-±S (12) 



S = Ucoth.a - i) (13) 

Denote by f{6) the inverse to 6. Since the 6 is odd and monotone function 
of (3, its inverse / exists and bears the same properties. Then the matrix B 
is the following function of the density matrix 



- — (^T',(4-))— /(I 
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Since the expression flT^ for p is the independence of the choice of 6, in two- 
dimensional case both matrices B and p are defined by their mean deviation 
values (3 and 5, respectively. So, the essential dependence of the matrix 
'temperature' parameter B from the density matrix p is completely captured 
by the function /. Its graph looks as follows. 




5 Lazy ensembles are equilibrium 

Like Gibbs ensembles in thermodynamics, the lazy ensembles are equilibrium, 
namely, the introduced 'temperature' parameters B possess the equalizing 
property. To show it, first introduce the notion of conditional ensemble. In 
terms of game played by Alice and Bob this means that Bob measures a fixed 
observable H upon the particles emitted by Alice. Again, he has the uniform 
distribution as null hypothesis, but the constraint in Q is of scalar rather 
than of matrix character. Solving the appropriate variational problem 



Sip) 



mm, 



P^pOj) dtp = TiHp 



we obtain 



(14) 



— this ensemble is conditional with respect to given observable H. 

Consider two quantum systems with state spaces Ti and 7i', respectively. 
Let their states initially be p and p' . Then, since we consider a non-interacting 
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coupling of the systems, the joint density matrix is p®p' in the tensor product 
space H ®H'. Let us measure the sum of values of the observables H and 
H', that is, introduce the observable H = i/®I'+I® H'. The conditional 
optimal ensemble of separable states with respect to the observable H is the 
following distribution 

M ^ /'^ [~/5h TrHP^55^/] 

i^M^-ip) = — 

Like in classical thermodynamics, the partition function of the joint system 
is the product of subsystems' partition functions: 




therefore the equalizing property holds 

If Ph < Ph' then Ph < Ph < Pw (15) 

which means that the conditional lazy ensembles are equilibrium and that P 
plays the role of temperature parameter. 

Concluding remarks 

Continuous ensembles of pure states proved their relevance in various as- 
pects of quantum mechanics. From the theoretical perspective, they provide 
the limit cases on which numerical characteristics of density matrices are at- 
tained, for instance, the minimal value of accessible information about the 
state is attained on 'Scrooge' ensemble which is a continuous distribution p. 
Furthermore, we claim that they are relevant from the operationalistic point 
of view. Even if we are speaking of preparing discrete ensembles, we must 
also have in mind that their are unavoidably smeared by various noises and, 
strictly speaking, we have to deal with continuous distributions. 

We use the techniques of continuous ensembles to carry out statistical 
inference in quantum realm according to the standard scheme: we have an a 
priori hypothesis (we necessarily need it, otherwise there is no way to make 
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any inference [7j), then we obtain some information about the system and 
have to shift to a new hypothesis. 

In our case the null hypothesis is the assumption that any pure state is 
emitted with equal probability. Then the information is obtained that the 
average density matrix of the state is p. We show how, starting from the 
'minimal effort' assumption, to guess the strategy of the preparation of the 
pure states. As a result, we obtain so-called 'lazy' ensembles. 

These ensembles are also proved to provide the minimal deviation from 
the null hypothesis. They are described by exponential distributions Q of 
pure states averaging to a given density matrix p: 

where the matrix parameter B plays a role in some respect similar to temper- 
ature, in particular, it is shown to possess the equalizing property. Although 
we may not treat it as a fully-fledged temperature, for instance, in contrast 
with classical thermodynamics, it is ambiguously defined up to an arbitrary 
additive constant. According to formula (jSj), we can so choose the additive 
gauge for B that In Z will vanish and the mean value Tr Bp will be equal to 
the differential entropy of the ensemble, so we may call this matrix parameter 
B 'differential entropy observable'. 
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